Collapse and revival of entanglement of two-qubit in superconducting quantum dot 
lattice with magnetic flux and inhomogeneous gate voltage 
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We study the entanglement of a two-qubit system in a superconducting quantum dot (SQD) 
lattice in the presence of magnetic flux and gate voltage inhomogeneity. We observe a universal 
feature for the half-integer magnetic flux quantum which completely washes out the entanglement 
of the system both at zero and finite temperature. We observe that the ground state is always 
in a maximally entangled Bell state when there is no inhomogeneity in gate voltage in the 
superconducting quantum dot lattice. We find an important constraint in magnetic flux for ground 
state entanglement. We also observe few behavior of entanglement at finite temperature is in 
contrast with the zero temperature behavior. 

PACS numbers: 03.65.. Ud, Entanglement and quantum nonlocality. 
PACS numbers: 74.78. Na, Mesoscopic and nanoscale. 
PACS numbers: 74.81. Fa, Superconducting wire networks. 



INTRODUCTION 



Recent years have witnessed a surge in research activ- 
ity involving a close interaction between the subjects of 
quantum information science and many body condensed 
matter physics [H-Q • The key concept in quantum 
information is that of entanglement, a truly unique 
feature of quantum mechanical systems. Entanglement 
implies non-local correlations between quantum particles 
which do not have a classical counterpart. Entanglement 
measures provide an additional characterization of 
the many-body state in condensed matter. Let us 
consider the state of an interacting spin system as an 
example. Since the ground state wave function changes 
qualitatively in a quantum phase transition, it is of 
significant interest to probe how the genuine quantum 
aspects of the wave function, namely entanglement, 
changes as the transition point is traversed. The 
question that arises is whether the entanglement of the 
quantum system extends over the macroscopic distances 
as ordinary correlation do. In recent years, quantum 
phase transitions have been extensively studied in spin 
systems using well known quantum information theoretic 
measures like entanglement, fidelity, reduced fidelity, 
fidelity susceptibility. The fidelity typically drops in an 
abrupt manner at a quantum critical point indicating a 
dramatic change in the nature of the ground state wave 
function 0-|l3. 

We have been motivated by the extensive studies of 
entanglement physics and interesting results in spin 
system [TlTj24T | . We decided to apply this concept of 
study in the different disciplines of quantum condensed 
matter quantum many-body system. Here we study 
the entanglement physics of two qubit superconducting 
quantum dot lattice. Before we proceed further, we 
would like to discuss the basic aspects of a superconduct- 



ing quantum dot lattice: Superconducting quantum dot 
lattice consists of array of superconducting grains. The 
superconducting grains are of nanoscale size, and differ- 
ent states of the superconducting grain are controlled 
by the ratio between charging energy and Josephson 
energy, while the average charge of the dot is controlled 
by the gate voltage. These are similar to the metallic 
and semiconductor quantum dot, where the different 
phases of the system depends on the ratio of Coulomb 
charging energy to the kinetic energy. Therefore the 
system can be described as a superconducting quantum 
dot (SQD) and the lattice as a superconducting quantum 
dot [25143 if lattice. 

One can consider the Cooper pair of SQD as a charged 
boson. The Physics of SQD can therefore be described 
in terms of interacting bosons. Bosonic physics is more 
interesting and hard to understand than the fermionic 
physics. One can understand the starting point of 
fermionic physics from the standpoint of independent 
electron approximation, at least in higher dimensions, 
whereas to understand the interacting bosonic system, 
one has to introduce the interaction from the very 
beginning. At the same time quantum phase diagram 
of this SQD lattice is very rich with different quantum 
phases. Therefore the study of the entanglement physics 
for the SQD lattice for both zero and finite temperature 
is interesting in its own right [25-27]. 
It is well known from our previous studies that 
mesoscopic SQUID array can also be treated as the 
superconducting quantum dot lattice with modulated 
Josephson junction [25-26]. The authors of Ref. [28| 
have found the magnetic flux induced superconducting 
Coulomb blocked in mesoscopic SQUID array and also 
the magnetic flux induced superconductor- insulator 
quantum phase transition. Experimentally and also 
theoretically it reveals that the applied magnetic flux has 



an important effect in the SQD lattice system. We will 
see in due course of our study that the inhomoginity of 
the gate voltage plays an important role in entanglement 
to disentanglement (product state) transition. Therefore 
we are motivated to study entanglement physics of SQD 
lattice in the presence of applied magnetic flux and 
inhomoginity of gate voltage. 

The plan of the manuscript is as follows. We present 
the model Hamiltonian and entanglement physics in 
section (II) of this manuscript. We present summary 
and conclusion in the section (III) of the manuscript. 



II. MODEL HAMILTONIAN FOR 
INHOMOGENEOUS SUPERCONDUCTING 
QUANTUM DOT LATTICE AND THE STUDY 
OF ENTANGLEMENT PHYSICS 

II A. Model Hamiltonian and Ground State 
Analysis 

At first we write down the model Hamiltonian of SQD 
lattice system with Josephson couplings having on-site 
charging energies and inter-site interactions in the pres- 
ence of gate voltage and external magnetic flux. We also 
consider the inhomogeneity in the applied gate voltage. 
The Hamiltonian is written as 

H = Hji + Heco + Heci- (1) 

We recast different parts of the Hamiltonian in quantum 
phase model as 

Hji = - Eji\cos(tv-§^)\ J2i cos(cf) l+1 - fa), 

where fa and fa+i are quantal phase of the SQD at the 
point i and i+1 respectively, as 

TT _ ECO V I i 9 2^2 

where Eco is the on-site charging energy, Now 
Heci — Ezi ^ rij n i+ i, 

i 

where Ezi is the NN charging energies between the dots 
respectively. In the phase representation, {—i-^p) is the 
operator representing the number of Cooper pairs at the 
ith dot, and thus it takes only the integer values (m). 
Here, Hamiltonian Heco accounts for the influence of 
gate voltage (eiV ~ V g ), where eN is the average dot 
charge induced by the gate voltage. When the ratio 
— ¥ 0, the SQD array is in the insulating state hav- 
ing a gap of the width ~ Ecoi since it costs an energy 
~ Eco to change the number of pairs at any dot. The ex- 
ceptions are the discrete points at N = (2n + 1), where 
a dot with charge 2ne and 2(n + l)e has the same energy 
because the gate charge compensates the charges of ex- 
tra Cooper pair in the dot. On this degeneracy point, a 
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FIG. 1: Color online, concurrence vs magnetic flux at zero 
temperature for different values of SV g = 0.0, 0.2, 0.5 for the 
dashed bullet and + signs respectively. Here Ej = 1.0, Eco = 
1.0 and Ezi = 0.4. 



small amount of Josephson coupling leads the system to 
the superconducting state. 

Here we recast our basic Hamiltonians in the spin lan- 
guage, where each site of the dot is either empty or singly 
occupied. During this process we follow Ref. [25[ and 
13 ■ Now 

H. n = -2 E J1 \co8{ir%)\'E i (SjS i+1 - + h.c), 

and 

Heco = ^£,(2^ -hf. 
Heco = — ^Eco^2nhiSi Z '. 

Here hi = jVi ~^"~ 1 allows the tuning of the system 
around the degeneracy point by means of gate voltage. 
We can tune the gate voltage in such a way that we 
can generate inhomogeneity in on-site charging energy. 
Without loss of generality we can also write the model 
Hamiltonian as 

Heco = ^(£co+<5^)^ z +Y,( E coSV g )S l+1 z . (2) 

i i 

Where 8V g is the variation of gate voltage around the 
lattice sites. 

Heci — 2Ezi J2i &i z Si+i z ■ 
The total Hamiltonian of the system is 

H = 2E J1 \cos(TT-^-)\Y,(SjS i+ r + h.c)+2E zl Y / S i 

+ ]T(£co + SVg)S* + E(£co - SV g )S i+1 z 

i i 

Now we consider the Hamiltonian for N = 2 case. We 
would like to write the Hamiltonian in the standard basis, 
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FIG. 2: Color online, concurrence vs. magnetic flux at zero 
temperature for different values of Josephson couplings. Here 
Ej = 1.0, SV g = 0.5 and E Z i = 0.5. 



the Boltzmann constant. We calculate the concurrence 
to measure the entangle of two qubit system of super- 
conducting quantum dot lattice following the reference 
of Wootters' formula [14,15]. The analytical expression 
for the concurrence is 

C = max (0, 2max\ l - A») (4) 

i 

, where A^ is the square roots of the eigenvalues of the 
matrix R = p{o\ v ®o 2 v )p*{o\ v ' ®o 2 v ). p* is the complex 
conjugate of p. The system is maximally entangled when 
C = 1 and the system is in disentangle state when C = 0. 
The density matrix of the system is 
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B = 2Eji\cos(tt-^)\. The eigenstates of this two 
sites Hamiltonian are >= |0, >, \ip% >= 

>H > = 
8V a -A 



|i,i >, 1^3 >= y5r*(ci|i,o > +lo, i >), |v> 

^H^i(pa|1.0 > +I * 1 >)• where ci 
A = 



4£ j V*(t£)I 



\{^Ezi — 2Eco), E 2 = ^(2Ezi + 2Eca), E 3 = —Ezi— A 
Ea = —Ezi+A. If we consider the homogeneous system, 
i.e., there is no variation of gate voltage over the lattice 
sites. The two states \tpa > and \ipi > are the maxi- 
mally entangled Bell states, i.e., (1/V2)(|0, 1 > — 11,0 >) 
, (1/V2)(|0,1 > +|1,0 >). As we see from our analyt- 
ical expression that ground state depends on the value 
of Eco i Ezi and A. Ground state is in the disentan- 
gle state (product state) when the ground state energy is 
either E\ or E 2 , otherwise the system is in the entangle 
state. Thus for this superconducting quantum dot lattice 
system there is a transition between the disentangle state 
to entangle state due to the variation of the system pa- 
rameters. We will see in due course of our study that the 
magnetic flux plays an important role in the transition 
between the disentangled state to the entangle state. 

II. B Entanglement Study for Zero and Finite 
Temperature 

Now we calculate the thermal entanglement of two ar- 
bitrary qubits (N = 2) of the superconducting quan- 
tum dot lattice. The density matrix of superconducting 
quantum dot lattice in equilibrium at temperature, T is 
p = ^exp(—H/ksT), where H is the Hamiltonian of the 
system, Z is the partion function of the system, ks is 



Ei = 



The square roots of the eigenvalues of the matrix R 
are Xi = \ 2 = \/A~[Dl : A3 = \/B 2 Cz + C 2 and A4 = 
y/B 2 C 3 - C 2 . Where A x = ±e- Ezi / T [cosh{E C0 /T) - 
^^sinhiEco/T) Di = ±e- E ^ T [cosh{E C o/T) + 

,zl- TS inh(Eco/T)} C 2 = -J^e^smhiA/T) 

B 2 = ±e Ez ^ T [cosh(A/T) - S -^-sinh(A/T)]; C 3 = 
±e Ezi / T [cosh(A/T) + ^sinh(A/T)}. In this derivation 
we consider ks — 1. 

Finally we obtain the formula for concurrence by using 
the relation Eq. 5, 

C = I e E ^ T (iEjainh^) e^l T ) (5) 

Now we calculate, the concurrence of the system by us- 
ing the above relation and the final expression for the 
concurrence at T = is 



C = \B/A\ 



(6) 



Now we analyze the above equation of concurrence in dif- 
ferent limit. If we consider the gate voltage in the lattice 
is homogeneous ,i.e., SV g — that implies the concur- 
rence C(T = 0) = 1. The system is always in maximally 
entangled state. For this situation magnetic flux has no 
effect to make a transition from entangled state to disen- 
tangled state. 

If we consider the presence of inhomogeneous gate volt- 
age, SVg 7^ then the system has finite concurrence but it 
has not maximally entangled for any arbitrary magnetic 
flux. We observe an universality that at half-integer val- 
ues of magnetic flux quantum, the concurrence is zero. 
The analytical expression which we derive for the con- 
currence is obtained satisfies the following relation 
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FIG. 3: Color online, concurrence vs. magnetic flux at differ- 
ent finite temperatures. Here Ej = 1.0, SV g = 0.5, Eco = 1-0 
and E z i — 0.5. 
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FIG. 4: Color online, concurrence vs. magnetic flux for dif- 
ferent inhomogeneity of gate voltage in SQD lattice. Here 
Ej = 1.0, E C0 = 1.0 and E zl = 0.5. 



Fig. 2 shows the variation of concurrence with magnetic 
flux for different values of Josephson couplings. We 
observe from our study that concurrence is large for the 
smaller values of ratio between the Josephson coupling 
to on-site Coulomb charging energies. As we have 
obtained from our previous studies that this large ratio 
of -J^ favors the Luttinger liquid and superconducting 
phase of the system and the smaller one favors the Mott 
insulating density wave phases of the system [25,26]. 
Thus it is clear from our study that entanglement of the 
system is large for Mott-insulating phase of the system 
where the Cooper pairs are localized. But the collapse of 
entanglement at the half-integer magnetic flux quantum 
has the universal feature for all values of Josephson 
couplings. 

Fig. 3 shows the variation of concurrence with magnetic 
flux at different temperatures. It reveals from our study 
that the concurrence is smaller for the higher values 
of temperature and at the same times the system is 
in the disentangled state for wider region of magnetic 
flux. The region of disentanglement state is larger for 
the higher values of temperature. The solid lines in the 
figures are for the considerations of superconducting 
Coulomb blocked effect induced co-tunneling effect in 
the system. We observe that co-tunneling effect has not 
drastic effect to change entangle transition of the system 
at the quantitative level. For ground state entanglement, 
the effect of co-tunneling effect is not explicit in the 
analytical expression (Eq. 8). 

Fig. 4 shows the variation of concurrence with magnetic 
flux for different values of SV g . Here we observe that the 
presence and absence of gate voltage inhomogeneity has 
no appreciable for the transition from entangle state to 
disentangle state. This behavior at finite temperature is 
in contrast with the behavior at zero temperature where 
the system is always entangled state for homogeneous 
gate voltage. 



From the analysis of the above equation we get the con- 
straint of magnetic flux that to obey the concurrence re- 
lation system. 

, ^ 4>o -i y/(E C o ~ 2E Z1 - SV g )^/{E C o - 2E Z1 + SV g ) 
a> > — cos — - 

7T 2Ej 

. (8) 

Now we present our results and physical explanations. 
Fig. 1, shows the variation of concurrence with magnetic 
flux for different values of 5V g . We observe from the 
figure that the system is maximally entangled Bell state 
when there is no inhomoginity in the applied voltage 
in the lattice. This maximally entangled Bell state 
is independent of the applied magnetic flux. We also 
observe from this figure that a finite inhomogeneity in 
the applied gate voltage introduces a transition from en- 
tangled state to disentangled state. In the entanglement 
phase system is almost maximally entangled state. 



III. SUMMARY AND CONCLUSIONS 

We have studied the entanglement of a two qubit sys- 
tem of a superconducting quantum dot lattice. We have 
observed the magnetic flux dependent entangled state 
to disentangled state transition at zero and finite tem- 
perature. We find an important constraint on magnetic 
flux during the study of ground state entanglement. We 
have observed that at zero temperature that the sys- 
tem is always maximally entangled state when there is 
no-inhomogeneity in the gate voltage in the lattice sites 
whereas the presence and absence of inhomoginity has 
no effect at finite temperature. We have found out the 
universal feature for half-integer magnetic flux quantum. 
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